Induced vortex tunneling in a superconducting wire 
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We consider induced topological transitions in a wire made from cylindrical superconducting film. 
During a transition, a pulse of electric current causes transport of a virtual vortex-antivortex pair 
around the cylinder. We consider both the instanton approach, in which the transition is viewed as 
motion of vortices in the Euclidean time, and the real-time dual formulation, in which vortices are 
described by a fundamental quantum field. The instanton approach is convenient to discuss effects 
of the environment, while in the dual formulation we show that there exists a potentially useful 
adiabatic regime, in which the probability to create a real vortex pair is exponentially suppressed, 
but the total transport of the vortex number can be of order one. 

PACS numbers; 85.25.Hv, 03.67.Lx 



I. INTRODUCTION 

Tunneling of magnetic flux across weak superconduct- 
ing links of various nature is of interest for both fun- 
damental physics and applications. In particular, de- 
vices using Josephson junctions have been discussed re- 
cently in connection with proposals for quantum comput- 
ing (qubits). (For a review of various types of supercon- 
ducting qubits, see ref. [1].) Such devices (SQUIDs) can 
be, for many purposes, reduced to just a few degrees of 
freedom, so it becomes possible to numerically simulate 
their quantum behavior. In this way, one can discuss 
not only the amplitude transfer between the basis states, 
but also the residual excitation due to population of the 
higher levels [2,3]. 

As the strength of the weak link is increased, by go- 
ing from a Josephson contact to a narrow wire, and then 
to a wider wire, superconductivity becomes more robust, 
and the tunneling rate drops. In ultrathin wires, with 
widths Lx ^, where ^ is the Ginzburg-Landau coher- 
ence length, the relevant tunneling events are quantum 
phase slips [4,5]. Crossover to the wide- wire case occurs 
at widths S,, and for wider wires transport of flux 

happens via formation of vortices. 

Robustness of superconductivity in wires with 
can be useful in applications that require long coherence 
time, such as qubits. Indeed, we expect that in a suit- 
able implementation such a wire will not be subject to 
any shunting effects except for thermal quasiparticles. 
Possible implementations include a wire defined litho- 
graphically on a substrate or a multiwall carbon nan- 
otube coated with a superconducting film. In either case, 
the superconducting circuit has to be closed to allow a 
persistent current, which presumably can be done using 
ordinary bulk superconductors. 

On the other hand, when a wire with Lj; 3> ^ is used 
in place of a weak link, inducing a transition between 
the basis states of the qubit becomes a nontrivial mat- 
ter, since flux does not travel easily across wider wires. 
We propose to use for this purpose a pulse of electric cur- 
rent along the wire, which can be produced by coupling 



the qubit, either directly or inductively, to some external 
circuit. The pulse will lower the potential barrier sepa- 
rating the basis states and thus encourage the motion of 
flux across the wire. We call this process induced vortex 
tunneling.* 

Because this process involves many degrees of freedom, 
the theory of it has to be constructed from an entirely 
different standpoint than the few-degree-of-freedom mod- 
eling common in the theory of SQUIDs. In this paper we 
develop and compare two approaches. One approach uses 
instantons, i.e., Euclidean solutions, which in our case 
correspond to motion of vortices and antivortices across 
the wire in Euclidean (imaginary) time. This approach 
makes it possible to discuss the effects of the environment 
(specifically, fermions at the vortex cores) on tunneling, 
but for calculations of the residual excitation left in the 
system by the pulse, is less convenient than the other ap- 
proach, based on duality. In the dual formulation, vor- 
tices are described by a fundamental quantum field, and 
the tunneling rate appears as a finite-size effect due to 
mode quantization in a confined geometry. 

For definiteness, we consider cylindrical geometry, in 
which the wire is formed by a thin superconducting film 
on an insulating cylindrical surface. The "width" of 
the wire is the circumference of the cylinder. Although 
this may not be the simplest geometry to manufacture, 
it has the advantage of simple (periodic) boundary con- 
ditions in the x direction. Results for other geometries, 
e.g., thin strips should be substantially similar. 

Short Abrikosov flux lines that can propagate in the 
film will be referred to as vortices and antivortices, de- 



*As we explain in more detail below, the pulse of current 
is supposed to enhance tunneling of fiux but not induce real- 
time, over-barrier transitions. Thus, our problem is quite 
distinct from that of a resistive state in thin wires [6] or films 
[7]. In particular, real-time classical simulations based on the 
time-dependent Ginzburg-Landau equation are clearly inap- 
plicable in our case. 
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pending on which direction the flux points. Motion of 
the flux across a cyUndrical wire can be viewed as creat- 
ing a vortex-antivortex pair on one side of the cyUnder, 
transporting them along the circumference, and annihi- 
lating them on the other side, see Fig. 1. This process 
has been discussed in various contexts in the literature 
[8-10] . Here, we present a detailed calculation of the rate, 
including the eSect of core fermions, and a discussion of 
the residual excitation. 




FIG. 1. A segment of wire made from cylindrical supercon- 
ducting film. A pulse of current I{t) induces transport of a 
virtual vortex-antivortex pair around the cylinder. 



apply duality at the second-quantized level. In a second- 
quantized description, vortices are represented by a fun- 
damental quantum field. Such a description was used in 
refs. [13,14] to study nucleation of vortex-antivortex pairs 
in the presence of a static supercurrent. The main differ- 
ences between that problem and ours are that in our case 
the process is, first, time-dependent and, second, avoids 
as much as possible production of real vortex pairs. 

The paper is organized as follows. Sect. II describes 
the duality map. Effects of core fermions and of the 
vortex-antivortex potential on vortex tunneling are dis- 
cussed, using the instanton picture, in Sect. III. This 
discussion gives us confidence that the linear equation (9) 
obtained in Sect. 2 indeed captures the essential physics 
of the process. Sect. IV contains definitions of the av- 
erage vortex current, which is used to measure the tun- 
neling amplitude, and of the residual excitation left in 
the system after the pulse (the number of real vortices 
produced). In Sect. V we consider the adiabatic limit 
(a slowly changing current) and show that it is possible 
to have a sizeable amplitude transfer with exponentially 
small residual excitation. Some numerical estimates are 
given in Sect. VI. We find that for a film of thickness 
d = 10 nm, a suitable width of the wire is of order 1 /xm. 
Sect. VII is a conclusion. 



Vortex instantons have been used to discuss the phe- 
nomenon of quantum creep [11]. That discussion em- 
ployed the ordinary non-relativistic vortex action. In 
our case, since a timneling event requires formation of a 
vortex-antivortex pair, we are constrained to use a "rela- 
tivistic" action, with the speed of light replaced by some 
other limiting speed ci. We will see that in this "rela- 
tivistic" case, instantons arc characterized by large (Eu- 
clidean) velocities and correspondingly short timescales, 
one of the consequences being that the core contribution 
to the vortex mass is greatly reduced. 

We stress that when wc speak about production of a 
vortex pair, transporting vortices around the cylinder, 
etc. we mean a virtual process, i.e., tunneling. Pro- 
duction of real vortices is detrimental to our goal, since 
these vortices would be easily "detected" by the environ- 
ment (e.g., by electrons at the vortex cores), and that 
would result in rapid decoherence of quantum superpo- 
sitions that we intend to form. To prevent decoherence, 
induced transitions should leave as little imprint on the 
environment as possible. The main result of the present 
work is that it is possible, at least theoretically, to use 
a slowly (a(liabatic;ally) changing current that has an ex- 
ponentially small probability to create a real vortex but 
still leads to a sizeable transport of the vortex number 
around the cylinder. 

At the first-quantized level, the dual description we 
use coincides with that employed in the theory of vortex 
transport in arrays of Josephson junctions [12]. How- 
ever, to consider effects of zero-point vortex-antivortex 
pairs fluctuating "in and out of existence", we need to 



II. DUALITY MAP 

The original quantum description of a superconductor 
contains the usual vector potential as well as fields 
corresponding to other degrees of freedom: fiuctuations 
of the order parameter, normal electrons at vortex cores, 
etc. We introduce field ax, dual to the electromagnetic 
field A^, by rewriting the corresponding factor in the 
functional integral as a Gaussian integral: 

exp I / F^„d^xdt \ = const, x 

I Va^ exp ^ J (^^e^'^^F.^a, + ga^) d'xdt . (1) 

Here d is the film's thickness. 



e is the magnitude of electron charge (e > 0), and c is the 
speed of light. We have already incorporated the condi- 
tion that the film is thin, so that parallel fields penetrate 
inside undiminished. Greek indices run over values 0,1,2; 
Latin indices over 1,2. Sums over repeated indices are 
implied. 

The next step is to separate the original field A into 
^(^rt)^ which consists of narrow flux tubes attached to 
vortices, and the remaining long-range part In 
thin films, the screening length for perpendicular mag- 
netic fields is very large [15], so a narrow flux tube does 
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not account completely for the magnetic field of a vor- 
tex: part of the field will appear in a^j]^ t^ih \)q 
responsible for long-range interactions between the vor- 
tices. However, the electric field of a moving vortex is 
confined to the vortex core [16,17,10], and that makes in- 
troduction of the flux tubes useful. For the field strength 
we have 



P — F(long) , p(vort) 

We define a (conserved) current of vortices as 

tA _ ff*" ^nv\ l?(vort) 



(4) 



Thus a state with Jq = C(5(x) carries precisely a unit flirx 
quantum $o = 27r/.g of F'^''°'^\ 

To obtain a dual description of vortices, we introduce 
a new complex field x with a unit charge with respect to 
the field a. The action of the vortex field is taken in the 
form 



-I 



(fxdt 



\dtX + {ic/h)aox\' 



-cl\diX+ii/n)aix\-M^{t)\x\' 



(5) 



where M{t) is half of the pair-production frequency. Dur- 
ing a pulse of current, this frequency is time-dependent. 
Eq. (5) is a "relativistic" action, except that the role of 
the speed of light is played by some other speed ci, cf. 
refs. [13,14]. We will discuss the value of ci in Sect. VI. 

The "relativistic" form of (5) may seem ad hoc, but in 
fact the form of the vortex action is highly constrained 
by symmetries. For example, the presence of a conserved 
vortex current dictates that x should be a complex field, 
while the possibility of pair production dictates that the 
kinetic term should be quadratic in the time derivative. 
The absence of terms with higher spatial derivatives is 
a good approximation as long as wavenumbers do not 
exceed a certain ultraviolet cutoff, which we expect to be 
of order of the lattice spacing. Finally, self-interaction of 
the vortex field (higher-order terms in x) can be neglected 
as long as vortices are sufficiently rare. 

The simple action (5) neglects other degrees of free- 
dom, in particular, fermions at the vortex core [18]. 
These are nearly gaplcss and, as well known, cause large 
friction for slow vortex motions [19]. For a tunneling 
process, involving virtual vortices only, the expected ef- 
fect of such a dissipative environment is a reduction in 
the tunneling rate [20] . We analyze this effect in the next 
section and find that for rapid, large-frequency Euclidean 
motions corresponding to tunneling, the role of dissipa- 
tive processes is greatly rcdiiccd, so much so that they 
do not affect the main exponential factor in the tunneling 
rate, but can only affect the preexponent. Because our 
main conclusions are based on the value of the exponen- 
tial factor, we take 



(6) 



as the full action of our dual description. 

The vortex current in this description is given by the 
usual formula = —6S-^/6a\. In particular, the spatial 
components are 



if 2i 



(3) From (6), we obtain an equation for a in the form 



(7) 



(8) 



Combined with eqs. (3) and (4) this is seen to be com- 
pletely equivalent to the saddle-point condition for the 
integral (1), i.e., to the original definition of a. 

Equation of motion for x is obtained from (5) and reads 

[dt + {ic/h)ao]\ - 4[di + {i/h)ai]\ + M\ = . (9) 



When we use (8) in this equation, we obtain an interac- 
tion of vortices with the long-range field F^°'^^\ (The 
interaction with the short-ranged is negligible if vor- 
tices are rare.) As we have mentioned above, the field 
includes the long-range tails [15] of vortex mag- 
netic fields and thus gives rise to a logarithmic inter- 
action between vortices. Accordingly, the frequency M 
in (9) should be understood as being due to the vortex 
core only. However, for instanton processes, involving a 
single vortex-antivortex pair, we can, instead of keeping 
track of the vortex contribution to a^, model the vortex- 
antivortex potential by making M position-dependent. 
Moreover, as we will see in the next section, for calcu- 
lation of the tunneling exponent, such x-dependent fre- 
quency can be replaced simply by its average value along 
the vortex path. Within this model, and with the short- 
ranged dropped from (8), components of a do not 
contain any fields due to vortices and reduce to external 
E and B fields only: 



(ao,ai,a2) 



2gc 



(Bz, -Ey, Ex) . 



(10) 



Eq. (9) is a second-quantized description of vortices, 
which allows us to study effects due to virtual vortex- 
antivortex pairs. Similar descriptions were used previ- 
ously in refs. [13,14] to study the instability of a static su- 
percurrent via production of real vortex-antivortex pairs. 
In the limit when vortex-antivortex fluctuations are ne- 
glected {not a suitable limit here), eq. (9) gives rise 
to a first-quantized description, in which a vortex of 
mass ra = HAl/cl moves under action of a Lorentz force 
built from dual electric and magnetic field. These fields, 
e = — Voo — (l/c)9ta and b = V x a can be expressed 
through the electric current and charge density, using 
(10) and Maxwell's equations. The resulting expressions 
coincide with those appearing in the first-quantized the- 
ory of ref. [12]. 
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For a static supcrcurrcnt / in the y direction, wc would 
have ai{t) cx It. In contrast, the adiabatic case consid- 
ered below corresponds to a slowly changing pulsed af. 
ai{t ±00) = 0. 

An additional constant "magnetic" field b can be used 
to describe a Magnus force on the vortex. However, as 
we will find in the next section, for tunneling processes 
the Magnus force is canceled, to a very high accuracy, 
by the spectral flow force [21-23] due to core fermions. 
Accordingly, we set the total effective Magnus force to 
zero: b = 0. 

Eq. (9) simplifies further if we restrict our attention 
to the cylindrical geometry of Fig. 1. A pulse of electric 
current will be applied lengthwise, in the y direction. 
The resulting field Ey will induce a virtual pair of vortex 
and antivortex to travel along the circumference of the 
cylinder (i.e., in the x direction). So, we set Ex = = 
0, and Ey = E{t). For the present configuration, E is 
independent of x and y. Thus, the equation of motion 
takes the form 



2 o2 



dx - ^-^E{t) 



X + M\t)x = 0. (11) 



Although wc will not need any such relation in what fol- 
lows, we note that the two time-dependent parameters 
in (11), E and M, can be expressed through a single 
function of time — the electric current (see ref. [24] for 
details). 



Situation, however, is entirely different at large fre- 
quencies, or short timescales, which as we will see below 
is the case relevant to tunneling. The characteristic fre- 
quency associated with the fermionic response is deter- 
mined by the "minigap" coq ~ A'^/ep, separating differ- 
ent states at the vortex core [18]. When the timescale of 
the vortex motion is much shorter than lOq^, the response 
and the contribution of fermions to the inertial mass will 
be diminished. Here, wc calculate this effect for the vor- 
tex motion corresponding to our tunneling problem. The 
motion takes place in the Euclidean (imaginary) time. 
We will employ the kinetic equation [23] for core states, 
which gives their evolution in real time, and then analyt- 
ically continue the results to the Euclidean domain. 

The low- lying ( "anomalous" ) branch of fermion spec- 
trum [18] at the vortex core runs in energy approximately 
from —A to A (for a review and further references, see 
ref. [26]). This branch is characterized by two quantum 
numbers: k^, the momentum in the z direction (which 
in our case is perpendicular to the film), and the angu- 
lar momentum I. For a vortex at rest, the energy lev- 
els are given approximately by Ei{kz) = —iOo{kz)l, with 
<^o{kz) = A'^ /2vFk±, and k"^ = kp — k^. In particular, 
Wo(0) = A"^ /2vFkF [23]. Motion of the vortex relative to 
the superfluid results in an interaction Hamiltonian and 
can cause transitions between different core levels. The 
distribution function n in general depends on the quan- 
tum numbers kz and I and their conjugate coordinates z 
and (j). If n is independent of z, the kinetic equation [23] 
reads 



III. INSTANTONS AND THE EFFECTS OF THE 
ENVIRONMENT 



dn 

'at 



oJo{kz)-g^+kx v{t) 



dn 
'dl 







(13) 



It is well known that motion of vortices in supercon- 
ductors is subject to strong renormalization effects due 
to the "environment" of normal electrons at the vortex 
cores. For low- frequency motions, in very clean super- 
conductors, the main effect is a drastic reduction of the 
limiting speed Ci in eq. (9). Indeed, in such cases it is 
natural to expect that c\ will be equal to the critical ve- 
locity of the superconductor Ucr ~ A/vp, where A is the 
value of the gap outside the core. At larger speeds, the 
vortex becomes a "tachyon" and causes an instability via 
production of quasiparticles from the core. Equivalently, 
the small ci = Vci can be interpreted as a large value of 
the inertial mass m of the vortex, since from (9)^ 



m ■■ 



M/c 



(12) 



This estimate is consistent with direct calculations [25] 
of the inertial mass (provided the frequency M is under- 
stood as the core contribution only). 



tin this section, we set h = \. 



where k = (fc_L cos^, fc_L sin*/), fcz), and v is the vortex 
velocity. Note that we write the kinetic equation in the 
collisionless (Vlasov) approximation. This is appropriate 
whenever the characteristic timescale of the motion is 
much shorter than the collision time. 

Solutions to eq. (13) can be obtained analytically for 
any form of the vortex velocity. Here we consider the case 
of main relevance to us, when the velocity is in the x di- 
rection at all times. Then, the zero-temperature solution 
is a step-function 



n{(p,l;t) 



l + kj_ujo / dt'GR{t-t')h{(p;t') 



where 



h{(l); t) = v{t) cos (t> + i^Q ^dtv{t) sin ( 



(14) 



(15) 



and Gn is the retarded Green function for df + ui^: 
Gn{t < t') = 0, Gnit > t') = uo^^ sin[wo(i - t')]. (More 
generally, the step-function can be replaced by an arbi- 
trary function in (14), so finite-temperature solutions can 
also be constructed.) Eq. (14) describes excitation of the 
core system, originally in its ground state, by the vortex 
motion. 
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If velocity v{t) is peaked around t = at the timescale 
ti shorter than both luq^ and the "observation time" t, 
we can replace it in (15) with a pulse of the form v{t) = 
x6{t), where x is the displacement. This yields n{<j), I; t) = 
0{l) for t < 0, and 



n((/), l;t) = 9[l + k±x sin(woi 



(16) 



for t > 0. Integrating over all modes, we obtain the 
total momentum transferred to electrons upon a vortex 
displacement x: 



Px{x) = 



J-kp 



(17) 



while for Py the sine gets replaced by cos Wot. The remain- 
ing dependence on t means that even after the displace- 
ment has been completed, and the vortex has stopped, 
fermions at the core continue to oscillate ("wobble"). 

Tunneling in a wire can be visualized as motion in Eu- 
clidean time of a vortex and an antivortex, each by half- 
width of the wire, or alternatively as motion of the vortex 
alone by the full width L^. For definiteness, we will dis- 
cuss the latter process. We refer to it as an instanton. 
The instanton has two timescales associated with it: time 
To it takes to create or destroy a vortex-antivortex pair, 
and time n it takes the vortex to travel the distance L^. 
Both To and ti are variational parameters that adjust 
themselves to achieve the largest tunneling amplitude. 

A natural estimate for tq is tq ~ M~^, and we will see 
that this estimate is unaffected by the presence of the 
environment. We will also find that in the limit when 
instantons are dilute, i.e., when the vortex motion has 
a large Euclidean action, we have n ^ M~^, but still 
wqTi <^ 1. So, to understand the effect of the environ- 
ment on Ti, we analytically continue to t = —ir and ex- 
pand in small loqt. The change in the vortex momentum 
is opposite to (17), and the corresponding contribution 
to the Euclidean action of the vortex is 



Pxdx = —LlkpU!oiO)T 



(18) 



In our case, we need to replace t with ti, and al- 
though this just invalidates the approximation under 
which (17) was obtained, and the precise form of the ve- 
locity pulse begins to matter, we can still obtain an order- 
of- magnitude estimate for J Pxdx. Let us first obtain this 
estimate for the case when the instanton gas just ceases 
being dilute, and ti becomes of order M~^. In the case 
of induced tunneling, parameters M and A should be 
taken at time when the transition actually occurs. How- 
ever, because (18) is now proportional to A^/M, and M 
scales as A^, we can as well use the unperturbed, equilib- 
rium values. Using M ^ (kFd)eF, and luo{0) = A^/4ei?, 
we estimate (18) as (fc_FLa;)^(A/ei?)-^/64. This is small 
for typical values kpL^ ~ 10^ and A/ep ^ 10~^, but 
can become large for larg er L/x- III the latter case, or in 



the dilute regime when ti is sufficiently large in compar- 
ison with M~^, the action (18) begins to play a role in 
determining the value of the variational parameter ti . 

To see how that happens, consider the Euclidean action 
for a single vortex, without the environment correction: 



JdrMil+viy/^ 



(19) 



where ve is the Euclidean velocity, and we use units 
with ci = 1. Strictly speaking, the mass M should be 
made cc-dependent, to include the vortex-antivortex po- 
tential. For now, however, we will consider it as some 
constant average mass; we will justify this replacement 
later. Then, during the vortex motion, ve is constant, 
and Ti = Lx/ve- The exponential factor in the tunneling 
amplitude in this case is given by 



dv 



Ve 



{l + vlf'^-\nvE 



(20) 



For large values of the product ML^, the saddle-point 
value of Ve is large: ve = {MLxY^^, while ri = 
(Lj/A/)^/^. The saddle-point value of the exponent in 
this case is 



Se = MLx/ci + 0(1) 



(21) 



(we have restored ci in this formula) , so the limit of large 
MLx is precisely the limit when the instanton gas is di- 
lute. Now, including the correction (18), with r ^ n, we 
obtain a new effective action, which in the limit of large 
Ve has the form 



5eflF 



MLx + ^^+CLlk%iJoiO)^+lnvE , (22) 



2vi 



16VE 



where C is a numerical coefficient of order 1. Depending 
on the parameters, either the second or the third term 
in (22) can be more important in determining the saddle 
point value of ve- In either case, however, the effective 
action on the saddle point will have the same form as 
eq. (21). We conclude that the longitudinal momentum 
transfer from the vortex to fermions at the core affects 
at most the tunneling preexponent, but not the main 
exponential factor. 

The exp(— const. La;) form of the tunneling exponent 
had been discussed in the literature for various systems 
supporting vortices [8,9]. The main point of this section 
is that it holds also for the "relativistic" action (19), and 
in this case the value of the exponent is not affected by 
the longitudinal correction (18). 

Turning to the transverse momentum, and following 
the steps that led to eq. (18), we obtain the following 
contribution to the Euclidean action: 



k^ d 

iPy{Lx)y = i-^LxV [l + 0(^0'^^)] 



(23) 
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where y is a smaU displacement in the y direction. Un- 
like (18), this contribution is purely imaginary: it gives 
rise to phase differences between vortex paths lying at 
different values of the lengthwise coordinate y. If signif- 
icant phase differences existed, the tunneling amplitude 
would be drastically reduced. However, the magnitude of 
the leading term in (23) is such that, for kp equal to the 
value of the Fermi momentum outside the vortex core, it 
precisely cancels the action corresponding to the Magnus 
force. This is the well known cancellation of the Magnus 
force by the force due to spectral flow [21-23,27,28]. The 
remaining correction in (23) is of order CiJq and is already 
quite small. Still, an even more complete cancellation 
can be achieved, because kp, which measures the elec- 
tron density at the core, is also a variational parameter 
and does not have to coincide with the value outside the 
core. As eq. (23) shows, a tiny self-adjustment of kp 
is sufHcient to render the efi'ective Magnus force exactly 
zero. 

Finally, we need to restore the effect of the vortcx- 
antivortex potential. As before, we consider the situation 
when, upon creation of a virtual vortex-antivortex pair, 
the antivortex stays in place, while the vortex makes a 
complete circle around the cylinder. Then, for the vortex 
action, we can still use eq. (19) but M should now be 
replaced by some (periodic) function of x. Accordingly, 
vp is no longer a constant. We notice, however, that the 
reason why (19) previously simplified into (21) was not 
that Vp was constant, but that it was large: vp ^ 1. 
This still applies now, during most of the path. Taking 
this limit in eq. (19) with M — > M{x), we obtain 



in the subsequent formulas, it is convenient to introduce 
instead of x a new field 



Then, 



X=i. (25) 



X^'da.X - d^X^'X - ^^E{t)X^X . (26) 



For our cylindrical configuration, the field X satisfies 
periodic boundary conditions in the x direction. The 
boundary conditions in the y direction will not matter; 
we take them to be periodic as well. The Fourier expan- 
sion of X is 



X(x,t) = {o^Mt) + Plj^) e*" 



(27) 



where k = {kx,ky), a and /3 are the usual annihilation 
operators, and the mode functions satisfy 



/k + o.^(t)/k = 
with a time-dependent frequency 



^i{t) = cik'y + ci 



k,--E{t) 



(28) 



+ M^{t) . (29) 



To achieve the correct commutation relation between 
X and dtX, the mode functions / should be normalized 
by the condition 



1 

Cl 



M{x)dx , 



(24) 



which is equivalent to (21) with M replaced by its average 
value along the path. 



IV. OBSERVABLES 

Wc have seen in the preceding section that the instan- 
ton picture of tunneling, with instantons corresponding 
to vortices and antivortices moving along closed paths 
in the Euclidean time, provides a very appealing and in- 
tuitive picture of flux tunneling across a wire of width 
Lx ^ ^. This picture has allowed us to estimate the 
effects of the environment and to obtain the main expo- 
nential factor in the tunneling amplitude. However, for 
certain tasks, in particular for understanding the adia- 
batic limit, it is more convenient to use the dual descrip- 
tion (9). In this section, we discuss what observables 
can be calculated in that dual theory that will contain 
information about tunneling. 

We will be interested primarily in two observables. The 
first is the average vortex ciirrent. The corresponding 
operator is given by (7); in our present case, only has 
a nontrivial average. To reduce the number of factors of h 



(30) 



where V is the total two-dimensional volume. At some 
initial time t = ti, before the pulse, we have E{ti) = 
and M{ti) = Mq, so we can use the usual plane- wave 
exponentials as initial conditions for f\^{t). Equivalently, 



/k(t^) = [2u;k(^^)^]-'/^ 
/k(ii) = -iuJkhiU) , 



(31) 
(32) 



where LoliU) = c^k'^ + M§. 

As discussed in the preceding sections, to include the 
effect of the vortex-antivortex potential, the frequency 
M{t) in these formulas is taken to be the average fre- 
quency along the vortex path. We note also that eq. 
(28) with the initial conditions (31)-(32) can be used for 
numerical studies of vortex transport [24] . 

Substituting (27) in (26), and assuming that the quan- 
tum state of the system is the vacuum of the operators 
a and /3, we obtain the average current as 



kx 



l/k 



(33) 



The integral of this over time gives the average total vor- 
tex number transported around the cylinder (i.e. in the x 
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direction), per unit length in the y direction. So, it mea- 
sures the efficiency of the amplitude transfer between the 
basis states of the qubit. 

The second quantity of interest is a measure of the 
residual excitation left in the system after the pulse. 
It is simply the total average number of vortices left 
at some final time t = tf and is obtained as the sum 
Sk"^k(i/) of the occupation numbers of all the individ- 
ual modes. These occupation numbers, as functions of 
time, are given by 



V 



\Mt)f + ujl{t)\Mt)\' 



(34) 



Our next goal will be to show that a sufficiently slow, adi- 
abatic change in E and M can lead to a sizable transport 
of the vortex number, while leaving nk{t ^ oo) exponen- 
tially suppressed. 



V. ADIABATIC LIMIT 

When the frequency (29) for each mode changes with 
time slowly (adiabatically), i.e. 



(35) 



the adiabatic theorem [29] guarantees that particle pro- 
duction will be absent in any finite order in l^tWkl/w^- In 
other words, nk(i), which was zero initially, will be zero 
at f ^ oo to exponential accuracy. The mode functions 
can then be approximated by WKB-type expressions: 



uj\^{t')dt' 



(36) 



Substituting this into (33), we obtain the average current 



as 



k 



where we have introduced notation 



E{t) ^ -E{t) . 



(37) 



(38) 



The summand in (37) depends on kx only in combina- 
tion kx — E, and it is odd in that combination. So, if we 
could replace the sum over kx by an integral and make 
a shift of the integration variable, we would prove that 
(37) is zero. There are, however, two obstructions to this 
procedure. First, the integral needs an ultraviolet regu- 
larization, which we take to be symmetric in kx, not in 
kx — E. As a result, the far ultraviolet modes contribute 
a finite amount proportional to E. Second, the difference 
between the sum and the integral results in a correction, 
which is periodic in E with period 2tt/Lx, where Lx is 
the circumference of the cylinder. This correction van- 
ishes when E is an integer multiple ofw/Lx, but is finite 



otherwise. It will be important for us to understand the 
structure of this correction. 

To make the argument more transparent, let us con- 
sider the case of a weak electric field, \E\ <C 2n/Lx, so 
that we can expand the expected periodic correction in 
E. To simplify things even further, we will also assume 
that 



, MLx 
6' = — ^ » 1 . 

2ci 



(39) 



Then, we can expand l/uiu in (37) in E, to obtain 
— {kx-E) 

Wk 

The first term in the braces gives zero upon summation 
over kx (assuming a symmetric ultraviolet cutoff), so we 
have 



iUt)) 



where 



cjE 
' V 



1 



{clkl + m2)3/2 



0(^2) 



m' = clkl + M'' . 



(41) 



(42) 



To compute the sum over kx in (41), we use the repre- 
sentation 



E 

n= — cso 

2 

in our case 



1 



(7r2n2 + 62)3/2 

1 



y/u)'^ — P sinh^w 



b = 



mLx 
2ci 



(43) 



(44) 



Under the condition (39), b is large, and we can expand 
the integral in (43) in e~^^. We obtain 



1 



^(C2fc2+ ^2)3/2 



1 + 2v^e-''' + 0{e-'^'') 



(45) 



If we were to neglect the discreteness of modes, i.e. 
replace the sum on the left-hand side of (45) with an in- 
tegral, we would obtain only the first term in the bracket. 
According to (41), the average current due to this term is 
proportional to E, with no other time-dependent factors 
(the factors of cancel out). In the absence of sources 
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of dissipation, such as pair-produced quasiparticles or 
real (non- virtual) vortices, we have E = —{l/c)dtA, so 
the integral of E over time is zero. Thus, the total vortex 
number transported around the cylinder is determined 
entirely by the exponential correction in (45). 

The remaining sum over ky can be replaced by an inte- 
gral and easily evaluated, noting that, at large values of 
6', e"^*" confines ky to rather small values: c\ky ~ M"^ /b' . 
Thus, in the limit (39) we finally obtain 



dtiMt)) 



1 



dtEMexp{-ML^/ci) . (46) 



This quantity has dimension of inverse length, as it gives 
the transported vortex number per unit length of the 
cylinder. 

The exponential factor in (46) coincides with the one 
obtained in sect. Ill by means of instanton calculus. 
Moreover, the periodicity with respect to E, noted af- 
ter eq. (38), has a simple explanation in the instanton 
picture. Indeed, from (11) we see that the product EL^ 
equals (up to a constant) the "solenoidal flux" f axdx. 
So, for ^ the amplitudes due to a single instanton 
and a single antiinstanton acquire opposite phases, equal 
in magnitude to \ELx\- In the limit (39), when the in- 
stanton gas is dilute, the average vortex current is then 
proportional to 



exp(— AfLj;/ci) sin EL^ 



(47) 



which has precisely the same periodicity as that deduced 
from eq. (37). We stress that in the dual calculation 
these results appear as a conseqiience of the quantiza- 
tion of vortex modes in a confined geometry, without 
any reference to instantons. So, the present calculation 
is complementary to the instanton calculation of sect. 
III. 

Note that oq. (46) determines also the value of the 
preexponent. As we have already discussed, unlike the 
exponent, the preexponent may in certain cases be af- 
fected by fcrmions at the vortex core. Any such modi- 
fication will be absent from eq. (46). We will, however, 
continue to use this equation, since our main conclusions 
arc based on the value of the exponential factor. 

The exponential suppression seen in (46) is of en- 
tirely different origin — and, hence, generally of different 
magnitude — than the adiabatic suppression of rik. To re- 
move the exponential suppression in (46) altogether, we 
need, at some time t during the pulse, to go just outside 
the limit (39), i.e. achieve the condition 



M{t) = ci/Lx . 



(48) 



(At this point, the 0{e^*'') terms in (45) will become 
important.) Eq. (48) will not jeopardize adiabaticity 
provided that the characteristic timescale tp of the pulse 
(e.g. the ramp time) satisfies 



On the other hand, one can choose parameters of the 
device so that, before and after the pulse, the proba- 
bility of random "errors", i.e., spontaneous transitions 
between the basis states, is vanishingly small. (The am- 
plitude of such spontaneous transitions is proportional 
to exp(— MqZ/x/ci), where Mq is the unperturbed initial 
mass.) 



VI. NUMERICAL ESTIMATES 

Let us summarize various conditions we have obtained 
so far on the parameters of the system. One is the in- 
equality (49), which is the condition that a significant 
(order 1) vortex number transport is compatible with the 
adiabatic suppression of production of real vortex pairs. 
Another is the condition of exponential suppression of 
spontaneous tunneling transitions ("errors"): 



» A = ci/Mo 



(50) 



The quantity A here, of the dimension of length, is anal- 
ogous to the Compton wavelength of an elementary par- 
ticle. 

To (49) and (50), we should add the condition that 
pair production of quasiparticles is also exponentially 
suppressed (since quasiparticles, like vortic;c!s. lead to 
decoherence-'-) and the condition that the coherence 
length of the superconductor remains smaller than L^, 
so that vortices remain well defined. These two condi- 
tions should be written so as to take into account the 
reduction of the superconducting gap A, and the corre- 
sponding increase in the coherence length ^, during the 
pulse: 



tp/2'K > h/A{t) 



(51) 
(52) 



Quantities without the bar denote the unperturbed val- 
ues, i.e. those before and after the pulse. We now obtain 
estimates for Mq, ci, and A. 

To obtain these estimates we use the Ginzburg-Landau 
(GL) energy of superconductor 



C|(V + i5A)V|^-|r||V|2 + ^|^|4 



(53) 



where r, and s are parameters. Although this ex- 
pression is strictly applicable only near critical temper- 
ature, we will use it at temperatures close to absolute 



ip/27r > 1/M{t) = L^/ci . 



(49) 



*A quasiparticle pair can detect passage of flux through 
Aharonov-Bohm scattering. Although it remains to calcu- 
late how strong this effect actually is, we nevertheless impose 
the condition of negligible pair production below. 
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zero, where our device is supposed to operate, since all 
we need are order-of-magnitude estimates. 

From the GL energy (53), Mq, which was defined as 
half of the threshold frequency required to produce a real 
vortex- antivortex pair, can be estimated as 



hMo ~ 27rdC|V'o| 



(54) 



Here is the order parameter in the absence of current. 
This is an estimate for tlie "core" frequency. The long- 
range interaction between a vortex and an antivortex will 
contribute an additional amount, enhanced relative to 
(54) by InLx/C However, since for the values of the 
parameters that we use below this logarithm is of order 
one, the estimate (54) will be sufhcient for our purposes. 
It can be rewritten using the London screening length 6: 



(55) 



where S^^ — 327re^C|'!/'oP/^^c^, and uem is the fine struc- 
ture constant. (In thin films, 5 determines the strength 
of the London current but not the actual screening length 
of the vortex magnetic field.) 

We now turn to estimating the limiting speed ci, or 
equivalently the vortex inertial mass (12). There are two 
types of contributions to m [27,17]: from the electric field 
caused by the vortex motion and from core fermions. We 
have discussed the mass due to core fermions in sect. Ill 
and have seen that for vortex tunneling, characterized by 
short timescales, the effect is not as large as for slow mo- 
tions: on a short timescale the vortex cannot transfer a 
large longitudinal momentum to the fermion subsystem. 
In particular, the core contribution to the mass does not 
affect the value of the instanton exponent. Accordingly, 
for our estimates we take m to be given by the small 
electromagnetic mass of the vortex 



m 



16e2^2 



(56) 



(See refs. [16,17,10] for how this estimate can be ob- 
tained.) This estimate results in a large value of ci: 



(57) 



The vortex's "Compton wavelength" can be now ob- 
tained from (50): 



d 



(58) 



For numerical estimates, we use the following values: 

d = 10 nm, ^ = 30 nm, and S = 100 nm. Then, according 
to (58), A ~ 35 nm, and (50) suggests that to suppress 
spontaneous transitions it is sufficient to use of order 
1 /im. Comparing (48) and (50), we see that in this case 
the pulse will need to reduce the vortex frequency M by 
a factor 



Mo 

M{t) ~ A ~ "^^ ■ 



(59) 



Because M{t) depends quadratically on the gap A{t) (cf. 
(54)), eq. (59) corresponds to a reduction in A, and an 
increase in ^, by a factor of order 5. For a superconduc- 
tor with a quasiparticle gap A = 10 K, this turns the 
condition (51) into 



tp/2iT > 4 X 10"^^ s 



(60) 



Meanwhile, the competing condition (49) is 

tp/27r > 10"^'' s . (61) 

So, for the above values of the parameters, adiabatic- 
ity with respect to quasiparticle production is a stronger 
condition than adiabaticity with respect to production of 
vortices. Finally, the condition (52) can be easily verified. 

VII. CONCLUSION 

Our main result is the demonstration that for topolog- 
ical transitions, involving transport of magnetic flux via 
vortices, there exists a potentially useful adiabatic limit. 
In that limit, the probability to produce real vortex or 
quasiparticle pairs is exponentially small (as might be ex- 
pected) , but the transport of the vortex number due to 
virtual vortices can in principle be of order one. 

Our choice of the cylindrical geometry was motivated 
by the ease of imposing bomidary conditions on the vor- 
tex field in the dual formulation. We expect the result to 
apply to other geometries, which may also be easier to 
fabricate. 

We have presented two calculations of the leading ex- 
ponential factor in the tunneling amplitude: one based 
on instantons, and the other based on the dual formu- 
lation, in which vortices are described by a fundamental 
quantum field. The results have been shown to agree. 

In addition, the instanton approach has allowed us to 
estimate the effects of the core fermions on vortex tunnel- 
ing. The process at hand, i.e. an adiabatic transition in 
which no vortices are produced in the final state, turns 
out to have a very short (Euclidean) timescale, result- 
ing in only a small transfer of longitudinal momentum 
to the fermions. As a consequence, the effective inertial 
mass of the vortex during tunneling is much smaller than 
the known estimates corresponding to slow motions. We 
have also discussed the effect of the vortex-antivortex po- 
tential and have found that, for calculating the leading 
exponential factor, the position-dependent vortex mass 
can be replaced with its average over the vortex path. 

Our calculations were done under the condition that 
the wire width (the circumference of the cylinder) 
exceeds the GL coherence length ^, so that vortices re- 
main well defined. (This is the limit opposite to that 
of an ultra-thin wire, where tunneling is due to quan- 
tum phase slips.) Our results suggest that for in the 
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micron range, induced vortex tunneling may be experi- 
mentally observable. 

Moreover, since the adiabatic limit, i.e. a slow switch- 
ing on and off of the current, eliminate the most obvious 
sources of decoherence, namely, pair production of vor- 
tices or quasiparticles, the process may be suitable for 
forming quantum superpositions of flux states. There 
are other possible sources of decoherence, in particular 
those due to components of the environment, such as 
nuclear or impurity spins, that are sensitive to magnetic 
fields caused by persistent currents, but we do not expect 
the situation here to be worse than for the conventional 
SQUID designs [30]. 

The author thanks Dane Bass and Albert Chang for 
discussions, and M. Stone and G. Volovik for correspon- 
dence on the physics of core fermions. 
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